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ABSTRACT 



We present an analysis of the large-scale molecular cloud structure and of the stability of clumpy structures in nearby molecular 
clouds. In our recent work, we identified a structural transition in molecular clouds by studying the probability distributions of gas 
column densities in them. In this paper, we further examine the nature of this transition. The transition takes place at the visual 
extinction of Ay' 1 = 2-4 mag, or equivalently, at E tal1 = 40 - 80 M Q pc~ 2 . The clumps identified above this limit have wide ranges of 
masses and sizes, but a remarkably constant mean volume density of 7? ss 10 3 cm 3 . This is 5 - 10 times larger than the density of the 
medium surrounding the clumps. By examining the stability of the clumps, we show that they are gravitationally unbound entities, and 
that the external pressure from the parental molecular cloud is a significant source of confining pressure for them. Then, the structural 
transition at Ay' 1 may be linked to a transition between this population and the surrounding medium. The star formation rates in the 
clouds correlate strongly with the total mass in the clumps, i.e, with the mass above A^' 1 , dropping abruptly below that threshold. 
These results imply that the formation of pressure confined clumps introduces a prerequisite for star formation. Furthermore, they 
give a physically motivated explanation for the recently reported relation between the star formation rates and the amount of dense 
material in molecular clouds. Likewise, they give rise to a natural threshold for star formation at Ay' 1 . 

Key words. ISM: clouds - ISM: structure - Stars: formation - dust, extinction - evolution 



1. Introduction 

Formation of dense, self-gravitating structures inside more dif- 
fuse, large-scale molecular clouds is the ultimate prerequisite for 
star formation. In addition to self-gravitating dense cores, molec- 
ular clouds in which star formation is taking place show ex- 
haustive structural complexity characterized by large contrasts in 
both density and velocity. From the general observation that al- 
most all known molecular clouds harbor young stars, it is known 
that the formation of structures capable of star formation (or 
alternatively, cloud dissipation) must proceed relatively rapidly 
compared to the life-times of molecular clouds. Likely as a re- 
sult of this complexity and rapid development, molecular clouds 
also show wide ranges of star-forming effic iencies and -rates 
(e.g. lHeiderman et al ll2010t lLada et"all2010h . This connection 
between the cloud structure and the capability of a cloud to form 
stars makes determining the roles of processes and parameters 
that control the cloud structure a fundament al open topic in the 
physics of star formation (review ed, e.g., by iMcKee & Ostrikerf 
l2007tlMac Low & Klessenll2004h . 

In the current analytic models of star formation, one par- 
ticularly important structural parameter of molecular clouds 
is the probability density function (PDF, hereafter) of vol- 
ume densities, which describes the probability of a vol- 
ume dV to have a density between [p,p + dp]. In such 
theories, the function has pivotal role: it is used to ex- 
plain among others the initial mass function of stars, and 
the star formation rates and efficiencies of molecular clouds 



(e.g. Padoan & Nordlundl l200l iKrumholz & McKee 2005: 



Elmeg reenll2008t Hennebelle & Chabrier 2009). In particular, 
this distribution is expected to take a log-normal shape in 
isothermal, turbulent media not significantly af f ected by the 
self-g r avity of gas (e.g. l Vazquez-Semadeni; ll994t IPadoan et ali 
Il997fc lOstriker et all Il999t iFederrath et al.l l2008bl) . Most im- 
portantly from the observational point-of-view, the log- 
normality of the distribution is expected to be reflected in 
the probabilit y distributions of column densities in molecu- 
lar clouds (IVazquez-Semadeni & G arcia 2001; Goodma n et al.l 
2009; IFederrath et al.ll2010l) . Also, recently a method has been 
developed to attain information of the actual three-dimensional 
density PDF b ased on the observ e d, two-dimen sional column 
density PDFs (iBrunt et alJl2010allbl : iBnmtl feOlO). Even though 
it has been pointed out that the general log-normal-like form 
for the (c olumn) density PD F can be borne out by various 
processes dTassis et al.l 1201 Oh . it is obvious that an account- 
able theory of cloud structure must meet with the observed 
characteristics of the distribution. This is particularly the case 
if the probability distribution shows any scale-dependent fea- 
tures and/or time evolution. Such properties have ind eed been 
predi c ted, e.g. in the presence of strong self-gravity dKlessen 
2000; IFederrath etaDboOSat ICho & Kind 120111: iKritsuk et al. 
2010), an d scale-dependent features have also been recently 
observed dKainulainen et alj|2009bt iFroebrich & Rowlesll2010l; 
iPineda et al.ll2010b ). This makes probing column density prob- 
ability distributions one measure of cloud structure that can be 
used to constrain analytic star formation theories. 
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However, the connection between theoretical and numerical 
predictions with observations of the PDF has been poorly inves- 
tigated. The studies in which the column density PDFs of mostly 
individual clouds have been examined have found a qualitative 
agreement with the predicted log-normal shape (e.g ., [Ridge et al. 
2006b} iGoodman et ail 120091; iBufler & Tan1l2009l) . The lack of 
systematic studies of column density PDFs has been mostly due 
to observational obstacles: all observational tracers of the cloud 
mass distribution suff er from shortcomings specific to the tracer 
in question (see, e.g.. IGoodman et aTl l2009). Generally, the dy- 
namical ranges probed by different molecular emission line trac- 
ers are often narrow, thereby probing only a limited range of 
the PDF. Dust continuum emission observations probe a wider 
dynamical range of column densities, but become insensitive at 
column densities below N < a few x 10 21 cirT 2 , thus missing 
a regime where most of the cloud mass is. This is the case also 
for dust extinction measurements using infrared shadowing fea- 
tures. In addition to these restrictions, mapping nearby cloud 
complexes that often span several degrees on the sky at high 
sensitivity requires a colossal observational effort, not generally 
feasible through typical observing campaigns. Dust extinction 
mapping in the near-infrared reaches only modest column den- 
sities of N < a few x 10 22 cirT 2 , thereby mostly missing dense 
star-forming clumps and cores. However, near-infrared extinc- 
tion mapping reaches very efficiently the low colum n densities 
N ~ 1 - 3 x 10 22 cirT 2 (Lombardi & Alvesl 12001). a regime 
where most of the cloud mass is. Therefore, it provides a feasi- 
ble tool to measure the column density PDFs at the scales of en- 
tire cloud complexes. The meth od has indeed been used rec ently 
for t his purpose, especially by [K ainulaine n et alj (l2009bl) (see 
also lLombardi et alJ 1200a l2008b| iFroebrich & Rowlesl l2010t 
Lomb ardi et alJl2010tlPmeda et al iGOlObl) ). 



In our recent work (Kai nulainen et al.ll2009bl Paper I here- 
after), we presented the first systematic study of the column den- 
sity PDFs in all nearby molecular clouds closer than 200 pc. We 
used near-infrared dust extinction maps of 23 molecular clouds 
to identify a transition in the PDF shape from a log-normal-like 
shape at lower column densities to a power-law-like shape at 
higher column densities. Such transition is characteristic to all 
star-forming molecular clouds. However, we showed that some 
of the non-star-forming clouds in our sample did not have the 
transition, but their PDFs were well fitted by a log-normal over 
the entire range of column densities above the detection limit. 
This led us to speculate that the PDF feature is linked to a transi- 
tion from a quiescent regime dominated by turbulent motions to 
a regime of active star formation dominated by gravity. The non- 
star-forming clouds not showing the feature would then belong 
entirely to the former regime, and during their subsequent evolu- 
tion towards star formation gravitationally dominated structures 
would appear, inducing also a transition to the shape of their 
PDFs. 

In this paper, we present a more detailed analysis of the struc- 
tural transition identified from the column density PDFs. In par- 
ticular, we will examine the physical characteristics and stability 
of the structures identified using the PDFs. With this analysis, we 
will show that the change in the PDF shape can be understood 
as a transition between the diffuse, interclump medium and a 
population of clumps that are gravitationally unbound, but sig- 
nificantly supported against dispersal by the external pressure 
imposed to them by the surrounding medium. This interpreta- 
tion links the PDF shape to a physically motivated explanation 
for the relation between star formation rates and the amount of 
high-density material in molecular clouds reported recently by 



lLada et al.1 (1201 Oh . and also to the threshold of star formation in 
molecular clouds. 

In §[2] we shortly describe the dust column density data used 
in this paper. In §f3]we characterize the structures identified from 
the column density maps and examine their stability. In §|4] we 
discuss the results and their impact for the structure- and star 
formation in molecular clouds . In {J5]we give our conclusions. 



2. The column density data of nearby clouds 

In Paper I, we use d the near-infrar ed color- excess mapping 
method presented b y Lombardi (120091) (see also lLombardil2 005 ; 
Lombardi & Alvesll200iir namely nicest, to derive dust extinc- 
tion maps for 23 nearby molecular clouds. The technique was 
used in conjunction with ne ar-infrared data from the 2MASS 
survey (Skruts kie et al.ll2006l) . resulting in dust column density 
maps covering the dynamical range of Ay ~ 1.2 — 25 mag 
in the spatial resolution of 0.1 pc (~ 2' at the distance of the 
nearby clouds). While the cloud sample for this paper is oth- 
erwise the same as in Paper I, we have excluded the Coalsack 
cloud from the analysis. Our recent molecular line observations 
of the Coalsack have shown that the region likely includes a sig- 
nificant extinction component not only from Coalsack, but also 
from an extended cloud at a larger distance (Beuther et al., in 
prep.). Since the effect of that component may well disturb the 
statistics derived in this paper, we decided to exclude Coalsack 
from the sample. 

As an example of our data, Figure Q] shows the extinction 
map derived for the Ophiuchus cloud. The figure shows also the 
PDF of the cloud, with the lognormal-like low-Ay part and the 
power-law-like tail at high-Ay clearly separable. The transition 
between these parts occurs approximately at A^' 1 ss 2.8 mag in 
this cloud. Throughout this paper, we refer to such position in 
the PDFs with A^' 1 . In the clouds included in the study, the tran- 

values, A*f = [2.0, 11] mag, 

tail 



sition occurs at relatively low Ay 

although in most cases between A^' 1 =2-4 mag. The A^' 1 values 
determined for each cloud are listed in Table [1] 

The A^ 11 value defines a set of spatially closed iso-contours 
in the column density maps (see Fig. [TJ. Throughout the paper, 
we will refer to the region where Ay < A^' 1 as the diffuse compo- 
nent, and similarly, to all regions where Ay > A'" 1 as the dense 
component. The former refers then, by definition, to the log- 
normal part of the PDF and the latter to the power-law-like part. 
Morphologically, the diffuse component is a uni-body structure 
in all complexes, but the dense component forms separate struc- 
tures. We will refer to such separate structures as clumps in this 
paper. 

3. Results 

In this section, we use the column density data introduced in 
§[2] to examine the nature of the diffuse and dense components. 
We first examine the physical characteristics of the components, 
namely the masses and sizes, densities, and velocity structure. 
Then, we consider the observed characteristics from the point of 
view of pressure balance in molecular clouds. 

3. 1 . Characteristics of the diffuse and dense components 
3.1.1. Total mass 

The gas column densities in the clouds can be inferred from the 
extinction maps using the measured extinction-to-gas column 
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Fig. 1. Dust extinction through the Ophiuchus molecular cloud, derived using the near-infrared extinction mapping technique and 
2MASS data (Paper I). The contours are drawn at A'" 1 = 2.8 mag, and at Ay = 1 mag which is used to calculate the total mass of 
the cloud. The mass in the dense component is ~ 21 % of the mass above Ay = 1 mag in this cloud (cf. Table [TJ. The inset shows 
the column density PDF of the cloud, with an approximate log-normal function fitted to the peak of it. 



ratio. We transformed the visual extinction values in the maps 
to hy drogen column densities using the relation dBohlin et al.l 
[19781) : 

/3 = (N(H) + 2N(H 2 ))/A V = 9.4 x 10 20 cm" 2 mag" 1 . (1) 

We then calculated the total mass of the cloud as a sum of ex- 
tinction values above Ay = 1 mag: 

Aftot = D 2 nP X J A v dxdy, (2) 

v7Q:Av>1 mag 

where D is the distance to the cloud, fi = 1. 37 is the mean 
molecular weight (adopting the same values as Lo mbardi et al.l 
(2008b), i.e., 63 % hydrogen, 36 % helium, and 1 % dust), and 
x and y refer to the map pixels. We adopt the same distances for 
clouds as listed in Paper I. We note that the chosen lower limit 
of Ay = 1 mag is arbitrary, and the total mass depends on the se- 
lected value (choosing lower threshold will yield higher masses 
for all clouds). However, a fixed value will make the values com- 
parable between the clouds. We also note that the Ay = 1 mag 
contour is closed in most mapped regions, thus uniformly defin- 
ing a cloud boundary. Table[T]lists the mean extinctions, Ay, for 
the clouds calculated using this definition for a cloud. 

The total mass of the dense component was calculated from 
the extinction in excess to the Ay threshold level: 

M dense = D 2 n/3x\ f A y dxdy-Afx f dxdy], (3) 

Jn:A v >A l f Jn:A v >A'f 



The mass of the diffuse component was then defined as Misuse = 
M tol - Mdense- We list in Table Q] the ratios of the mass of dense 
component to the total mass of the cloud. Clearly in all clouds, 
the mass of the diffuse component dominates the cloud mass. 
The ratios vary from a few percents for clouds with low star- 
forming activity to ~20 % for the most active clouds. Note that 
the total mass of the cloud was calculated as the mass above 
Ay > 1 mag. Since the column density below this level is, 
of course, not zero, our total masses represent lower limits. 
Accordingly, the quoted Mfe Dse /M tot ratios represent upper lim- 
its. 

3.1 .2. Clumps in the dense component 

In the following, we characterize the individual structures, i.e. 
clumps, in the dense component. We used a simple threshold- 
ing approach to identify the clumps from the extinction m aps, 
namely the "clumpfind2d" routine dWilliams et al.l [19941) . All 
pixels in the map that are connected with each other and above 
Ay 1 ' are considered as one clump. We emphasize that we do not 
make an effort to identify single-peaked structures nested in- 
side the contour defined by A^' 1 , because we particularly want 
to examine the mean physical parameters inside regions defined 
by the A^' 1 threshold. Therefore, there can be numerous distinct 
column density peaks (of any column density higher than Ay' 1 ) 
nested inside the clumps. In terms of the clumpfind2d algorithm, 
this approach equals to using A^' 1 as a threshold level for struc- 
ture detection, but not defining any additional column density 
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Table 1. Molecular clouds and the derived properties. 



Cloud 




4 90a 




Mot [10 4 M G ]* 




^clumps 


Star-forming 


clouds, physical resolution 0. 1 pc 






Ophiuchus 
Taurus 


2.8 
4 


6.7 
7.7 


2.4 
2.2 


0.52 
0.99 


0.21 
0.09 


20 
33 


Serpens^ 
Chal 


2.1 


3.2 


1.7 


0.27 


0.16 


11 


Chall 


2.1 


3.2 


1.8 


0.11 


0.17 


7 


Lupus III 


3 


4.1 


1.7 


0.08 


0.07 


18 


CrA cloud 


2.3 


4.1 


1.7 


0.10 


0.13 


11 


Lupus I 


2.0 


3.3 


1.6 


0.13 


0.12 


16 


LDN1228^ 


2.1 


3.1 


1.5 


0.15 


0.09 


4 


Pipe c 
LDN134 rf 


2.2 


3.7 


1.5 


0.09 


0.07 


3 


LDN204 d 


3.2 


5.3 


1.7 


0.26 


0.04 


10 


LDN1333 rf 


3.1 


5.0 


1.5 


0.25 


0.01 


9 


Non- star-forming clouds, physical resolution 0. 1 pc 


LDN1719" 






1.8 


0.18 






Musca 


2.0 


3.2 


1.7 


0.04 


0.16 


3 


Cha III 


2.1 


5.5 


1.7 


0.12 


0.11 


26 


Lupus V 






1.8 


0.28 






Star-forming clouds, physical resolution 0.6 pc 


Ori A CMC 


3 


5.6 


2.2 


9.2 


0.20 


15 


Per cloud 


3 


4.9 


2.0 


1.4 


0.12 


10 


Ori B GMC 


2.1 


5.4 


1.9 


7.0 


0.19 


35 


Cepheus A 


3.6 


13 


2.7 


1.1 


0.10 


3 


California 


4.2 




1.7 


10 


0.01 


14 



clouds in agreement with Larson's relations (M oc R 2 ). Larson's 
mean-density size relationship, p oc R' 1 , however, has been 
questioned by Ballestero s-Paredes & Mac Lowl (f2002). suggest- 
ing that it is an observational artifact due to limited dynamical 
range in column density. Our mass-radius relation, M oc r 2 - 7±0 - 2 
can also be seen as an indicator of the fractal dimension of 
the cloud, D « 2.7 ± 0.2. This range is consistent with the 
lar gest values in the range, D = 2.3 - 2.7, previously found 
bv lElmegreen & Falgarondtf 19961). indicating fa i rly space-filling 
colum n density structures (ISanchez et al.ll2005t Federrath et al.l 
I2009I 

The mean volume densities in the clumps, defined as n — 
is on the order of 10 cm . The mean volume 
densities of all clumps in all clouds have a distribution that peaks 
strongly at n « 0.8 x 10 3 cm 4 (shown in Fig. 0. The peak of 
the distribution is also relatively narrow, with the mean density 
being between n = 0.4 - 2.1 x 10 3 crrT 3 for 90 % of the clumps. 
We also calculated the mean density of the diffuse component 
in the clouds. This is straightforwardly defined for each cloud 
without the clumpfinding process, since the diffuse component 
is always a uniform structure. Similarly with the density of the 
clumps, the volume was calculated using the effective radius R = 
y/A/n for the cloud. The resulting mean densities, listed also in 
Table Q] are «diff = 1 — 2 X 10 2 cirT 3 , i.e. 5-10 times smaller 
than the mean densities of the clumps identified from the dense 
component. 



" Extinction above which the contribution of the log-normal compo- 
nent to the PDF is less than 10 %. 
* Calculated using the column density values above Ay > 1 mag. 
c could not be defined properly for the cloud. 
d The most prominent Lynds Dark Nebula in the region. 

levels that would be used in detecting peaks inside this parental 
structure. 

This exercise resulted in identification of ~ 10 clumps from 
each cloud complex that show wide ranges of sizes and masses. 
The effective radii of the clumps, defined as R = -\JA/tt where 
A is the area, varies between ~ 0.1 - 3 pc. Figure [2] shows the 
size distribution of all clumps in all clouds, showing that smaller 
regions are always more numerous than larger ones down to 
the resolution limit of our data {R as 0.1 pc). The size distri- 
bution has a power-law-like shape with the approximate slope 
of -0.9 + 0.2. Figure [2] also shows the mass distribution of the 
clumps, calculated by integrating Eq. [3] over the clump. The 
mass distribution covers roughly four orders of magnitudes be- 
tween 10 _I - 10 3 M Q , with approximately a power-law dis- 
tribution that has a slope of -0.4 + 0.2. This slope is flatter 
than typically observed f or the mass distributions of cores in 

the clouds ( 1 .3, e.g. iMotte et alj [19981: lAlves et alj I200H 

I Andre et"aTll2010h . being closer to the slopes derived for molec- 
ular clouds_or_CO c^himps_wijlrm_j^ 

e.g. IWilliams et aOl995l: iKramer et al.lfl998l: iBlitz et alj[2007h . 
It is, however, possible that the derived slope is affected by the 
blending of clumps with each other. Such blending can make 
the detection efficiency of clumps a function of the clump mass, 
and thereby affect the slope of the observed mass function (e.g., 
Kainulaine net~ai]l2009al: iPineda et alJ [2009). The mass-to-size 
relation for the clumps, also shown in Fig. [2] follows approx- 
imately the relation M oc r 2j±0 - 2 ^ again being cl ose to wha t 
has been derived for CO clumps within clouds (e.g.. lCarri ri987). 
The relation is also in agreement with what is expected for con- 
stant mean volume density spheres and steeper than expected for 



3.1.3. Correlation with CO and linewidth 

In Fig. [3j we demonstrate how the spatial extent of the dense 
component compares to the common molecular line tracer ob- 
servations. We use as an example the Ophiuchus and Perseus 
clouds for which large-scale 12 CO and 13 CO data are publicl y 
available through the COMPLETE survey (Rid ge et al1 l2"b06a). 
Figure [3] shows the 12 CO total antenna temperature map of the 
Ophiuchus cloud with a contour of Ay 11 = 2.8 mag overplotted. 
In Ophiuchus, thresholding at Ay l] separates two larger clumps: 
the main cluster region, and the streamers leading east from the 
cluster (there are additional clumps identified outside the cover- 
age of the CO emission). In general, the A^' 1 contour coincides 
quite well with the extent of the 12 CO line emission data (1-cr 
rms error of CO data is 0.98 K), while 13 CO is spatially less ex- 
tended. Figure [3] also shows a similar comparison on a smaller 
spatial scale for the B5 globule in the Perseus cloud. 

We identified 10 clumps in the Ophiuchus and Perseus 
clouds that are fully within the region covered by the 
COMPLETE survey. We estimated the virial parameters of these 
clumps, defin ed as the ratio of kinetic -to-gravitational energies 
in the clump (iBertoldi & McKedll992h : 



where G = 1/232 M ' pc (km s~') 2 is the gravitational con- 
stant and cr the velocity dispersion. The linewidths were esti- 
mated from both the 12 CO and I3 CO data by calculating the 
mean spectrum over the clump and making a simple gaussian 
fit to the peak of it. The mass was calculated from the extinc- 
tion data following Eq. [3] This calculation yielded virial param- 
eters a = 3 — 100 for the clumps. The virial parameters correlate 
with the mass of the clumps approximately in a power-law fash- 
ion (Fig. 0J. A simple linear least-squares fit to the data points 
yields the slopes of -0.69 + 0.12 and -0.64 + 0.13 for 12 CO and 
I3 CO, respectively. This relation is consistent with the prediction 



4 



J. Kainulainen et al.: From chaos to confinement 




-1.5 -1.0 -0.5 0.0 0.5 1 

Log R [pc] 



1 000.0 ^ c ) 
100.0 
10.0 



1«R~ 



•• t I +- ++ 

1 -o k kfe#+ 



00 



0.' 



.0 



I0.0 



R [pc] 



1 r 



N°<M~ 




-1 



2 3 



Log M [M ] 




000 2000 3000 4000 5000 



Fig. 2. Characteristics of the structures (clumps) defined by a thresholding at A^' 1 (see 93.1 .ll in text), a) The size distribution. Error 
bars show the y[n uncertainty. The dotted line shows the least-squares fit to the distribution, with the slope -0.9 ± 0.2. b) The same 
for the mass distribution. The linear fit results in the slope -0.4 ± 0.2. c) The mass-radius relation. Overplotted are slopes indicating 
M oc R} and M oc R 2 (dotted lines), the linear fit to the data points which has a slope 2.7 ± 0.2 (dashed line), and the resolution limit 
R = 0.1 pc of the data (solid line), d) The distribution of mean volume densities. 



for clumps con fined by ambient pressure from the medium sur- 
rounding them (Bertold i & McKeell 19921) and such relation has 
been previously observed for clumps identified from CO emis- 
sion data (e.g.. IBertoldi & McKedll992t I Williams et alJll995t 
lLadaetal.ll2008l) . It is. however, clear that the determination 
of the virial parameters suffers from likely non-gaussian errors, 
arising most pressingly from the uncertainty in determining the 
linewidth that would well trace most of the gaseous material in 
the cloud. Therefore, we consider this observed correlation in- 
dicative, although clearly not well constrained. It is, however, 
evident that having a » 1, these clumps are not gravitationally 
bound entities (although they can be significantly supported by 
other forces, as will be discussed later). This is unsurprising, as 
CO clumps in molecular cl ouds are generally observed to have 
high virial parameters (e.g., ICarrll987i IBertoldi & McKeell992t 
Falg arone etaD 19921) . 

Figure |4] also shows the size-linewidth relation for the 
clumps in Ophiuchus and Perseus. The data are scattered with 
no clear correlation, although similarly with the virial parameter 
- mass relation determining the correlation is hampered by the 
small number of clumps. We note that if the size and linewidth 
are uncorrected, the relation between the virial parameter and 
mass established above (a oc M~ 2 ^) implies that M oc R 3 , i.e. 
the mean volume densities of the clumps are constant (see Eq. 
@). Indeed, this is in agreement with the characteristics of the 
clumps derived in 93.1.21 i.e. that the mass-radius relation fol- 
lows approximately M oc R 2 - 7±0 - 1 , and that the mean volume den- 
sities of the clumps strongly peak around a characteristic value 
of - 10 3 cm 4 . 



3.2. Pressure confinement of molecular clumps 

As demonstrated in 93.1.11 the mass in the clumps above Ay 1 , 
i.e. in the dense component, accounts only for the minor frac- 
tion of the total gaseous mass of a cloud. In other words, the 
clumps are surrounded by a medium whose total mass (and spa- 
tial extent) greatly exceeds that of their own. As an example, 
the most massive clump in our cloud sample has the mass of 
about 10 % of the mass of the whole cloud (for other clumps, 
the fraction is much smaller). Likewise, as shown in 93.1.21 the 
mean density of the clumps is close to an order of magnitude 
higher than the mean density of the medium surrounding them. 
Thus, it seems reasonable to consider the gravitational force of 
the surrounding medium as a source of external pressure sup- 
po rting the clumps. In the fo llowing, we follow the formulation 
of IBertoldi & McKed (119921) and examine the scale of external 
support provided to the clumps by the diffuse medium surround- 
ing them. 

The basic condition for the virial balance of a clump is: 



= W + 2(7 - T ext ) + M. 
In this equation, W is the potential energy: 
3 GM 2 



(5) 



(6) 



and T and r ext are the kinetic energy of the clump and its surface 
term: 



T = 3/2 x P kin V, 



(7) 
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Fig. 3. Comparison of the dense component, i.e. structures above Ay > Ay 11 , with the CO molecular line emission. Top: 12 CO 
line emission from the COMPLETE survey for Ophiuchus, with a contour of A^ 1 ' = 2.8 mag overplotted. Bottom row: Similar 
comparison for the B5 globule in Perseus. The left panel shows the extinction map, with black contours at Ay = [1,3] mag. The 
whit e contour shows the extent of the coherent core in which the linewidth of the NH3 molecule emission drops abruptly, identified 
bv lPineda et alj (|2010a). The center and right panels show the same for the 12 CO and 13 CO line emission, respectively. 



r ext = 3/2 x p exl v. (8) 

M is the magnetic energy which we neglect for simplicity. The 
virial balance equation expressed in terms of pressure is then: 

Pan = Pp + Pext- (9) 

In this, the total kinetic pressure of the clump, i\; n , is the sum of 
both thermal and non-thermal components: 

Pkin = P(Ot + CT 2 T ), (10) 

and P gl is the gravitational pressure of the clump supporting it 
against expansion: 

Pgr = -1/3 xW/V = (^)G(pP) 2 . (11) 

Pext is the pressure external to the clump. Under the assump- 
tion that molecular cloud complexes are close to gravitational 



virial equipartition (e.g.,|Larso nll98lUHever et al.ll200ll) . a sup- 
porting external pressure is directed to a clump. This pressure 
arises from the turbulent pressure that balances the cloud against 
its own gravity. Since the cloud, as a whole, is close to virial 
equipartition, the turbulent pressure amounts to the gravitational 
pressure of the cloud (analogously to Eq. ITTt . but we adopt 
a slightly modified e xpression that takes into account that the 
cloud is not spherical dBertoldi & McKed 19921) . With the defini- 
tion of the mean mass surface density, S = M/(nR 2 ), the external 
pressure supporting clumps against dispersal is: 

Pext = Pfr d = (^p)G2 2 G , (12) 

where a\ and are numerical constants related to cloud 
morphology whose value can be evaluated as prescribed in 
Bert oldi & McKed d!992h . As an example of the order-of- 
magnitude of these pressures, using the typical 13 CO linewidth 
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Fig. 4. Top: Virial parameters derived for clumps identified in 
the Ophiuchus and Perseus clouds using thresholding at Ay = 
Ay 11 . In calculating the virial parameter, linewidths derived from 
13 CO data from the COMPLETE survey (iRidge et al.ll2006ah 
were used. The dashed line shows the slope corresponding to 
the relation a oc M~ 2 ^, predicted for c lumps confined by exter- 
nal pressure (Berto ldi & McKeel [l992). A linear fit to the data 
points yields the slope -0.64 ± 0.13. The open circles show the 
modified virial parameters of the same clumps, a, that take into 
account the external confining pressure (see 83.21 in text). The 
dotted line indicates the a — 1 (and a — 1) level. Bottom: Size- 
linewidth relation for the same clumps, shown for both 12 CO 
(red) and 13 CO (blue). Th e dotted line shows the cr oc R 05 rela- 
tion dSolomon et al.lfl987h . not a fit to the data points. 



of cr = 0.75 km s for a R = 1 pc sized clump, n^iff = 150 cm 
"clump = 800 cm 4 , and the mean mass surface density A 



-3 



v = 2 

mag yields the pressure ratios P^ n m 10 x P gr * 4 x f e xt- In 
other words, the pressures supporting the clumps against disper- 
sal amount in total to about one third of the pressure driving their 
dispersal. 

In the following, we examine these pressures for a popula- 
tion of clumps whose properties equal to those derived for the 
clumps identified in this paper. In 83.1.31 we showed that the ob- 
served virial parameters of the clumps scale with their masses 



(which is predicted for clumps confined by external pressure, 
iBertoldi & McKeeff9 92): 



2T 

W\ 



rjdump 

kin 
rjdump 



C\ X M 



-2/3 



(13) 



where c\ is a proportionality constant. It directly follows from 
this dependence that the internal kinetic pressure of the clumps is 
only a function of their density (Eqs.[T3land[lT1i. and thereby the 
kinetic pressure is constant for a population of constant density 
clumps. The ratio of outwards to inwards pressures for a clump 
is then: 

->clump 

(14) 



rjdump 



^ext + Pv 



(p) 2 ' 3 R 2 



9ai0G (2) 2 



16 (p) 4 ' 3 

Figure [5^ illustrates this ratio (Eq. [T4l as a function of the 
mean mass surface density X (in units of Ay) for clumps of dif- 
ferent sizes (R = 0.1-2.1 pc). As the mean density, we used the 
value n = 800 cirT 3 shown earlier to be the peak of the mean 
densities in the clumps ( 83.1.2b . Figure [5] shows that the transi- 
tion from a regime where structures are unbound to a regime 
where they are bound occurs around Ay ~ 4 mag (S a 80 
M Q pc~ 2 ) for clumps over a wide range of sizes (R = 0.1 — 2.1 
pc). This value is larger than the observed mean mass surface 
densities by a factor of about ~ 2 (see Table Q]), and therefore, 
the external pressure is lower than the internal kinetic energy of 
the clump by a factor of ~ 4. 

We have so far assumed that the kinetic energies of the 
clumps scale with mass as shown by Eq. [13] In principle, this 
scaling is supported by the linewidth data of clumps showing 
rather constant linewidths (Fig. [4}. This observation is, however, 
hampered by the poor statistics we could achieve with the avail- 
able data. Therefore, we consider also a case where the kinetic 
energies of the cl umps scale according to a Larson-like size- 
linewidth relation (Solom on et al.l [T987): 



.clump 



2T 

W\ 



pclump 5 x Q j2(Rf -5R 



kin 



rjdump 



GM 



(15) 



Again, Fig. shows the ratio of outwards to inwards pressures 
for clumps (Eq. fT4l as a function of mean mass surface density 
of the cloud. Using this scaling of kinetic energies, the balance 
occurs approximately at the level of the observed mean mass 
surface density (Ay ~ 2 mag). While our observations seem to 
favor kinetic energy scaling with mass (see also discussion in 
§|4]i , this example illustrates the behavior of the pressure balance 
in another plausible scaling scheme, suggesting that the external 
pressure indeed can be close to the internal kinetic energy of the 
clumps. 

Finally, we illustrate the net effect of the pressures to the in- 
dividual observed clumps by calculating modified virial param- 
eters for the clumps for which we have CO data. In this modified 
virial parameter, we take into account the pressure external to 
the clumps: 

p clump 



— clump 



kin 



, clump 



(16) 



The modified virial parameters are shown in Fig.|4]together with 
the traditional virial parameters that include only the gravita- 
tional and kinetic energies of the clumps. In agreement with the 
earlier results, the modified virial parameters are clearly smaller 
compared to the ones resulting from Eq. |4] but still somewhat 
larger than unity. 



7 



J. Kainulainen et al.: From chaos to confinement 




0246802468 
A (cloud) [mag] A (cloud) [mag] 



Fig. 5. a) The ratio of the pressures supporting a clump against a collapse (total kinetic pressure) to the pressures promoting it 
(external and gravitational pressure) in clumps with n = 800 cirT 3 . The different curves are for clumps with radii R = 0.1 - 2.1 pc, 
in steps of 0.2 pc. The dashed vertical lines indicate the interval of the observed mean extinctions (see Table [T). b) The same, but 
using Larson's size-linewidth relation to calculate the kinetic energies of the clumps (i.e., Eq.[l5]instead of Eq.[T3T>. 



4. Discussion 

4.1. Pressure confinement of the clumps 

We described in the previous section a new approach to char- 
acterize structures observed in molecular clouds. In particular, 
we used the observed gas column density PDFs of molecular 
clouds to define a population of clumpy structures (dense com- 
ponent) embedded in the extended, interclump medium (diffuse 
component). The transition between these components occurs 
at the extinction threshold A^' 1 =2-4 mag, or equivalently, at 
2 tail =4O-8OM pc~ 2 = 0.008 -0.017 gem 2 . This level is rel- 
atively constant, being in the quoted range in every cloud except 
one (Serpens, cf. Table [TJ for which it could be reliably deter- 
mined. The dense component becomes dominant at Ay = 3-8 
mag. The mass of the dense component is between 1 - 20 % of 
the total mass of the cloud which we defined as the total mass 
above Ay > 1 mag. The clumps of the dense component show 
roughly power-law-like distributions of sizes and masses, cov- 
ering wide dynamical ranges (see Fig. However, the com- 
ponents are characterized by remarkably constant mean volume 
densities of n ~ 10 3 cm -3 and n m 1.5 x 10 2 cm -3 for the dense 
and diffuse components, respectively. 

The clumps identified using the column density PDFs of 
the clouds in this study are very similar to the gravitational ly 
unbo und 13 CO clumps identified in seve r al studies in the past 
(e.g. ICarill987t iBertoldi & McKedfT99l IWilliams et aflll995t 
lLada et all 20081) . In particular, the mean densities, and mass- 
radius and virial parameter-mass relations derived in 93.1.21 
13.1.31 are in agreement to what has been derived for such CO 
clumps. Similarly, those studies have concluded that the external 
pressure can be a significant confining source for such clumps. 
A simple qualitative comparison between the extent of the re- 
gions selected with A^' 1 and 13 CO emission indeed suggests 
these methods may trace quite similar components in the clouds 
(see Fig. [3}. 

It has b een earlier suggested (iGoodman et al.l Il998t 
ICaselli et al.l 12002). and more recently directly observed 
dPineda et alj|2010al) that there appears to be a sharp transition 
to dynamically coherent objects, or cores, at the scale where 
non-thermal motions cascade from the supersonic to subsonic 
regime (the sonic sc ale, e.g. IVazquez-Semadeni et al.1 [ 2003; 
iFederrath etall 12010b . In particular. iPineda et al.l (1201 Oal) de- 
tected such transition in the B5 globule in the Perseus cloud, 
at the length scale on the order of R * 0. 1 pc, in agreement with 
the sonic scale (and Larson's size-linewidth relation). The spatial 



extent of this coherent core is illustrated in Fig. [3] together with 
the extent of the clump defined by Ay 11 . The transition to coher- 
ence occurs clearly in a different column density regime than the 
break in the column density PDF defined by Ay 1 '. Importantly, 
while the structural transition from supersonic- to subsonic ve- 
locities in cloud structure seems to be linked to a particular size- 
scale, the structural transition described by the Ay' 1 threshold 
is not size-dependent. As shown in Fig. [2] structures above A^' 1 
cover a large size- and mass range, with their number decay- 
ing in a roughly power-law-like fashion. Thus, the A^ 1 ' thresh- 
old appears to be unrelated to the transition to coherent cores in 
the velocity structure. Unfortunately, the column density maps 
used in our work do not provide high enough spatial resolution 
to properly sample the PDF at the length scale of the transition 
to coherence. Therefore, we could not directly look for possible 
features the transition would induce to the PDFs. 

We showed in Section 33.21 that a significant external pres- 
sure from the surrounding cloud is imposed to the clumps 
we identified using the A^ 1 ' threshold. This is especially the 
case if, instead of a constant kinetic pressure, clumps follow a 
Larson-like size-linewidth scaling relation. This indicates that 
the clumps may even be close to a pressure balance with their 
surroundings. The CO linewidths and virial parameters we de- 
rived for a small sample of clumps partially support this pic- 
ture: virial parameters (as defined by Eq. |4) correlate with 
clump masses as predicted for pressure confined clumps, and 
the linewidths do not show clear correlation with clump sizes 
(we observe a nearly constant linewidth for the clumps). On the 
other hand, the modified virial parameters taking the external 
pressure into account (Fig.|4]and Eq.[To*ll were somewhat in ex- 
cess to unity for those clumps, implying that they may be over- 
pressurized and thus either an additional pressure component 
may be significantly affecting them, or they may be expanding. 
This result is , again, sim ilar to what has been derived for 13 CO 
clumps (e.g.. lCarrlll987l) . 

The role of the internal gravitational pressure of the clumps 
is further illustrated in Fig. [6] which shows the mean density 
for all identified clumps as a function of the clump size (plus 
signs in the figure). A constant ratio of gravitational-to-external 
pressures defines a linear relationship in this plot with a slope of 
-1 (cf. Eq.fTTIandfLTt. The relation corresponding to P gI = P ext 
is overplotted in Fig. [6] All clumps identified from the dense 
component have densities lower than this relation, implying that 
the gravitational energy is indeed small compared to the external 
pressure. The typical mean density n = 150 cirT 3 of the diffuse 
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component is also shown, which obviously is clearly below the 
mean densities of the clumps. 

Despite the limits imposed by the spatial resolution and dy- 
namical range of the extinction maps, we can still examine the 
role of gravity in the structures nested inside the clumps (i.e, in 
smaller-scale structures inside what we have defined as a clump). 
We show in Fig. [6] with red diamonds a population of structures 
identified by an experiment in which we defined clumps using a 
threshold level A v = 3 x A 1 ™ 1 . While this threshold is typically 
Ay = 6-12 mag, such selection likely represents a popula- 
tion connected to star-forming regions, or at least, pre-stellar ob- 
jects. The structures identified with this experiment are above the 
Pgr = Pext lin^B This demonstrates how gravitation becomes an 
increasingly important confining force for density enhancements 
nested inside the clumps. To illustrate one case where gravita- 
tion is known to eventually become the dominant force, we have 
marked into the diagram the star-forming clump B5 in Perseus 
(see Fig. [3}. The clump defined by thresholding at A'" 1 is marked 
with black filled circle, and the structure identified inside it with 
the threshold at 3 x A^' 1 is marked with red filled circle. 

Given these results, we suggest that the observed organi- 
zation of structures, identified with a new approach using the 
column density PDF, can be understood as a population of 
clumps significantly supported by the external pressure. This ex- 
ternal pressure originates from the turbulent pressure outside the 
clumps, and in this framework, it is a consequence from the as- 
sumption that clouds as a whole are close to virial equipartition. 
Then, the break observed in the column density PDFs at A^ 1 ' rep- 
resents a transition from a diffuse inter-clump medium to clumps 
significantly supported by external pressure. This interpretation 
has some profound implications. Most pressingly, it implies that 
the external pressure from the large-scale cloud has an important 
role in the formation of molecular cloud structures over wide 
si ze and mass scales . This result is analogous to the recent work 
of lLada et al.l d2008l) . who found the external pressure a signifi- 
cant force in confining small-scale cores, or globules, in the Pipe 
Nebula. We note that in Paper I, we speculated that the transition 
in the PDFs could be related to a transition from a regime where 
gravitation is negligible to a regime where it is dominant. The 
results of this paper, however, revise this picture and more quan- 
titatively connect the break to the pressure conditions of density 
enhancements with their surroundings. 

We emphasize three main assumptions of this framework. 
First, we assume that molecular clouds are in virial equilibrium 
as a whole, in order to relate their gravitational pressure to the 
turbulent pressure confining the clumps. Although close to be- 
ing in virial equilibrium, molecular clouds are also embedded 
in an external medium and may have been formed in large-scale 
converging flows, being highly dynami cal objects, where surface 
terms play an important role (see e.g.,lBallesteros-Paredesl 2006: 
Bane riee etaT1l2009t ISchneideret al.ll2010l) . Second, as can be 
seen from Fig. [5] the external pressure due to the turbulent pres- 
sure has a significant role in providing support for the clumps. 
The confining force exerted to a clump by it, however, depends 
on the isotropy of the turbulent flow surrounding the clump. 
Supersonic turbulence is highly anisotropic locally and will lead 
to transient formation and destruction of filaments and clumps. 
While estimating the net effect of the turbulent pressure in a 



1 We note that this relation is not in any dependence to the pressure 
balance or virial status of the clumps (see Eg. II 1 land I12t. In order to 
estimate the kinetic pressure of the clumps identified in this experiment, 
a tracer probing densities of these objects would have to be used instead 
of 12 CO or 13 CO. 



more consistent way for each clump would require more detailed 
velocity information (and modeling of the cloud structure) than 
are available for this work, we estimated the average level the 
pressure may have on clumps implicitly in Eq.Q~2] Third, we ne- 
glected the magnetic field stress in our approach. The exact role 
of the magnetic fields in shaping cloud morphology is still un- 
der debate, but generally, it can provide confining pressure per- 
pendicular to the field lines which supports clumps against both 
collapse and expansion. The magnitude of this support can be 
estimated roughly with B 2 1 {(In), which generally leads to pres- 
sures simila r to turbulent ra m pressures for the field strengths of 
B * 15 fiG dCrutcherll999l). 

Interestingly, Crutc her et al.l (1201 Ol) recently found that at 
densities lower than n < 300 crrT 3 the magnetic field strength 
does not scale with density, implying that below that density the 
cloud material is likely channelled along the field lines. Above 
this threshold density, the field strength approximately scales 
with density as B <x n° 65 . In the context of our work, the thresh- 
old density of n w 300 cirT 3 is in the regime between the diffuse 
component (H w 150 cm 4 ) and the clumps (n » 800 cm 4 ). 
This rises the interesting possibility that the break in the column 
density PDFs at A^' 1 would be related to the change in the role 
of the magnetic field from dominant at low colu mn densities to 
less significant at higher column densities. While ICrutcher et all 
d2oioir suggests that this could be the regime where cloud struc- 
tures become self-gravitating, our work rather suggests that this 
change in the B-n relation could be due to a transition to struc- 
tures that are not self-gravitating but confined due to their ap- 
proximate pressure balance with their surroundings. 

4.2. Pressure confinement and star formation 

When coupled with the main result of Paper I, i.e. that non- 
star-forming clouds do not exhibit similarly strong tails (if any) 
in their PDFs as all star-forming clouds do, the interpretation 
discussed in this paper leads to a picture in which the forma- 
tion of pressure confined clumps occurs in clouds prior to (or at 
clearly higher rate than) the formation of gravitationally domi- 
nated cores. Indeed, in our sample of molecular clouds, pressure 
confined clumps are observed in some clouds that do not show 
active star formation (Musca, Cha III), or even high column den- 
sity cores (Musca, Hacar et al. in prep.). This picture is further 
supported by the recent analysis of the stabil ity of dense core s 
in the nearby, mostly quiescent Pipe Nebula (lLada et al.l l2008). 
In this cloud, Lada et al. examined a sample of ~ 150 cores of 
masses between 0.2 - 20 M Q {R w 0.04 - 0.2 pc). They con- 
cluded the core population in the Pipe to be pressure confined, 
gravitationally unbound entities. The cores in Lada et al. study 
were defined to be single-peaked (or at most a-few-peaked) en- 
tities, and thereby they may well represent the smallest scale of 
the hierarchy whose largest scale is represented by the structures 
identified in our study are. A similar result highlighting the role 
of extern al pressure in a star -forming cloud was recently pub- 
lished by Mar uta et al.1 ([2010). They investigated the stability of 
dense cores in the Ophiuchus cluster, in the regime of R < 0.1, 
and concluded that the external pressure has a significant role in 
the dynamics of the cores. These results clearly indicate that the 
hierarchy of structures nested inside clumps is affected by the 
external pressure all the way to the regime of dense star-forming 
cores. 

In the interpretation discussed above the formation of pres- 
sure bound clumps can be seen as a prerequisite for the forma- 
tion of gravitationally bound cores. Given this, there evidently 
should be a relation between the occurrence of such clumps 
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Fig. 6. Relation between the mean density and size of the struc- 
tures (clumps) defined by thresholding at Ay 11 . The dotted line 
shows the typical density of the interclump medium (i.e. the den- 
sity of the diffuse component). The dashed line shows the den- 
sity above which the gravitational energy of the clump becomes 
larger than the external surface pressure (P gl = P ext ). The red di- 
amonds show the structures identified from the column density 
maps by thresholding at Ay = 3 x A^ 1 '. The resolution limit of 
the column density data is at log R = - 1 .0. The black solid circle 
marks the clump B5 in Perseus, and similarly, the red solid circle 
marks the same clump when using the threshold of Ay = 3 x A^' 1 . 



and star formation, even beyond the general observation that not 
all quiescent clouds show such clumps, while all star-forming 
clouds do. Therefore, it is interesting to consider the observed 
star-forming efficiencies and -r a tes in the clouds of our sample. 

Recently, iHeiderman et al.l (l2010h studied the star-forming 
activities of nearby molecular clouds as a function of the gas 
surface density (i.e, the Kennicutt-Schmidt law). In their work, 
Heiderman et al. used near-infrared extinction maps similar to 
those employed in this paper to derive gas surface densities. 
They examined the number of young stellar objects (YSOs) in 
the clouds identified using the Spitzer satellite data in different 
column density intervals and constructed the Kennicutt-Schmidt 
law for their cloud sample. In particular, they observed an abrupt 
drop in the star formation rate at £ » 5 - 100 M pc~ 2 (see 
Figs. 3 and 8 in IHeiderman et afll2010l) . leading them to sug- 
gest a threshold for star formation at 2 t h = 129 ± 14 M pc~ 2 
(Ay _ 8.6 mag). A very similar result was reached recently by 
lLadaetalJd2010l) . who examined the relation between the num- 
ber of YSOs in nearby clouds and the amount of high column 
density material in them. They showed that the correlation be- 
tween the mass of the gas and the number of YSOs identified 
in the clouds is strongest (i.e., the dispersion in the relation is 
smallest) at A K ~ 0.8 mag (A v > 7.3 mag, or E » 1 16 M pc -1 ). 
We note t hat the dispersion of the SFR-surface density relation 
derived bv lLada et"aU (|20 1 Ol) starts to decrease already at surface 
densities lower than 2*116 M pc -1 , reaching its minimum at 
that point. 

The star formation thresholds derived in the studies above 
are slightly larger than the typical A^' 1 values. However, we de- 
fined the A^ 11 value as the point where the dense component, on 
average, becomes a significant excess over the diffuse compo- 
nent. Obviously, at such surface density the largest contribution 
to the PDF still comes from the underlying diffuse component, 
not from the excessive dense component. Typical surface den- 



sity values at which the contribution of the tail to the PDF be- 
comes dominant (> 90 %) are around A V J1 (90%) ~ 3 - 8 mag 
(listed in Table Q]). Such values would be v ery much in agree- 
me nt with the threshold v alues derived by lLada et al. (2010) 
and lHeiderman et ail d2010l) . given the very different approaches 
used in these papers. Therefore, it seems plausible to interpret in- 
crease in star-forming activity to be related to the regime where 
the column density PDF is becoming completely dominated by 
the dense component. 

In the context of clumps bound by external pressure, a natu- 
ral threshold for star formation is introduced by the surface den- 
sity at which pressure bound clumps form, which is around Ay 11 . 
Furthermore as discussed above, the mass above this threshold 
is in a direct connection to the SFR of the cloud, with the SFR 
increasing in a power-law manner with increasing gas surface 
density. This interpretation gives a physically motivated expla- 
nation for the star formation threshold occurring at relatively 
low surface densities and links it to an observed structural fea- 
ture in the clouds. Thus, we suggest a picture in which the forma- 
tion of pressure bound clumps, and thereby the structural transi- 
tion at Ay , introduces a prerequisite for star formation, with the 
amount of mass in clouds above that limit directly proportional 
to the capability of the cloud to form stars. 



5. Conclusions 

In this paper, we presented an analysis of the large-scale, clumpy 
structures in nearby molecular clouds and of their stability. In 
particular, we described a new approach to identify structure in 
clouds using the observed column density PDFs. With this ap- 
proach, we identified two distinctive components in them, re- 
ferred to as the dense and diffuse components, and described 
their basic physical characteristics. We then examined the sta- 
bility of the clumps in the dense component, especially by con- 
sidering the scale of external pressure imposed to them by the 
medium surrounding them. The main conclusions of our work 
are as follows: 

1 . The transition between the diffuse and dense components oc- 
curs at a narrow range of column densities, Ay" 1 = 2-4 mag, 
or S tal1 = 40- 80 M pc~ 2 . The dense component domi- 
nates the observed column density PDFs above Ay > 3 - 8 
mag. The total mass of the dense component is 1-20 % of 
the total mass of the cloud, and thus always clearly smaller 
than the mass of the diffuse component. Clumps identified in 
the dense component show wide dynamical ranges of sizes 
(0.1 - 3 pc) and masses (10 _1 - 10 3 M ). However, the 
mean volume density of the clumps is remarkably constant, 
n w 10 3 crrT 3 . This is ~ 5 - 10 times larger than the mean 
volume density of the diffuse component, n m 1-2x1 2 
cirT 3 . 

2. The clumps identified using the column density PDFs are 
gravitationally unbound and the external pressure, caused 
by the turbulent pressure from the diffuse, large-scale cloud 
surrounding them, can provide significant support for them 
against dispersal. However, examination of the stability of a 
small sub-sample of clumps indicates that they may be over- 
pressurized and either expanding or additionally supported 
by a component not included in our analysis (e.g. magnetic 
field support). Then, the physical properties of the clumps 
resemble those of the clumps often identified from 13 CO 
emission observations as structures of the lowest hierarchical 
level. 
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3. In lKainulainen et alJ (f2009b). we showed that some non-star- 
forming clouds do not show the PDF break, while some of 
them show a weak break but no gravitationally dominated 
dense cores. Coupling those results with the physical charac- 
teristics of the clumps derived in this paper suggests a picture 
in which pressure confined clumps form prior to, or at higher 
rate compared to, the formation of gravitationally dominated 
dense cores in the clouds. This suggests that the formation of 
pressure confined clumps is a prerequisite for star formation, 
and introduces a natural threshold for star formation at A^' 1 . 

4. The star formation rate in the cloud complexes of our sam- 
ple correlates strongly with the mass in the structu res defined 
by the A^ 1 ' threshold, as pointed out recently by lLada et all 
(2010), and furthermore, drop s abruptly below that surface 
density dHeiderman et alJfeOlOl) . This supports the interpre- 
tation laid out in Item 3 above. Most importantly, the inter- 
pretation then provides a physically motivated explanation 
for the relation between star formation rate and the amount 
of den se m aterial in the cloud s reported by Heiderma n et al.l 
d20loh and lLada et alJ d20loh . 
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